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ABSTRACT 

We have re-analyzed archival Ginga data of the soft X-ray transient source 
GS 2023+308 (V404 Cyg) covering the decline phase of its May 1989 outburst. Our 
spectral modeling includes the relativistically smeared Compton reflected continuum 
and iron Ka fluorescent line near 6.5 keV produced by X-ray illumination of the 
accretion disc. This gives a powerful diagnostic of the accretion geometry, with the 
amplitude of the reprocessed components showing the solid angle subtended by the 
disc, while the detailed shape of the relativistic smearing shows how close this material 
is to the black hole event horizon. 

We find that reflection is always significantly present in the spectra, but that its 
fractional contribution decreases as the decline progresses. The amount of relativistic 
smearing is also consistent with decreasing during the decline, although this is poorly 
constrained except for the brightest spectra. One plausible scenario explaining this 
evolution is of an optically thick disc with inner radius increasing as a function of 
time, with the X-ray source in the form of a central corona. This is similar to the 
evolution inferred for other X-ray transient sources, such as Nova Muscae, except 
that the underlying power law spectrum of GS 2023+338 stayed constant as the disc 
geometry changed. This challenges the underlying assumption of almost all models 
for the spectra of accreting black holes, namely that the hard X-rays are formed by 
Comptonization of seed photons from the accretion disc. 

Key words: accretion, accretion disc - black holes - stars: individual (GS 2023+338) 
- X-ray: stars 



1 INTRODUCTION 

The 'no hair' theorem for black holes states that they can 
be described solely in terms of their mass and spin (and 
charge, although this is unlikely to be important in phys- 
ically realistic situations). Black holes are then completely 
characterised by two parameters, making them much less 
complex than e.g. neutron stars, where there can be a strong 
magnetic field, and where there is an uncertain equation of 
state for the dense material. Black hole binary systems are 
then the simplest objects in which to study the physics of 
accretion flows in strong gravity, and have the further ad- 
vantage that the orbital parameters are often well studied 
so that the inclination, mass and distance of the system are 
tightly constrained. 

* Present address: NASA/Goddard Space Flight Center, Green- 
belt, MD 20771, USA 



Many of the black hole binary systems show dramatic 
outbursts where the luminosity rises rapidly from a very 
faint quiescent state to one which is close to the Edding- 
ton limit, Z/Edd, and then declines approximately exponen- 
tially over a period of months. These systems are the soft 
X-ray transients (hereafter SXT; see Tanaka & Lewin 1995, 
Tanaka & Shibazaki 1996 for reviews), with the outburst 
triggered by the classic disc instability which occurs as Hy- 
drogen starts to become partially ionised (Meyer & Meyer- 
Hofmeister 1981; Smak 1982; see Osaki 1996 for review). The 
dramatic increase in mass accretion rate onto the central 
object gives rise to strong X-ray emission. This illuminates 
the disc, and controls its evolution during the decline phase 
(King & Ritter 1998, King 1998), and in quiescence (King, 
Kolb & Szuszkiewicz 1997). 

The SXT also show spectral evolution during the out- 
burst and decline. In general, when the luminosity is close 
to the Eddington limit, their X-ray spectra show a strong 



© 1999 RAS 



2 P.T. Zycki, C. Done and D.A. Smith 



thermal component at ~ 1 keV, accompanied by a strongly 
variable power law (very high or flare state). As the lumi- 
nosity decreases then the power law component decreases in 
importance, giving the high or soft state which is dominated 
by the thermal component. There is then a marked transi- 
tion at L ~ 0.11/Edd as the thermal component decreases 
in both luminosity and temperature, and a hard, strongly 
variable power law dominates the X-ray spectrum (the low 
or hard state). This spectrum remains fairly stable to very 
low luminosities, where the source is in the quiescent or off 
state (Miyamoto et al. 1992; Nowak 1995; Tanaka & Lewin 
1995; Tanaka & Shibazaki 1996). 

According to current paradigm the thermal component 
originates in an optically thick accretion disc for which the 
Shakura-Sunyaev (1973; hereafter SS) solution is commonly 
used, while the hard power law is produced by hot, opti- 
cally thin plasma. The geometry inferred for the low/hard 
state from modelling the broad band (1-200 keV) spectrum 
of both transient and persistent sources (Dove et al. 1997; 
Gierlihski et al. 1997a; Poutanen, Krolik & Ryde 1997; see 
Poutanen 1998 for a review) is that of a geometrically thin, 
cool accretion disc which truncates at some transition ra- 
dius, r tr , which is rather larger than the innermost stable 
orbit at 6R g (R g = GM/c 2 ), with the inner region filled by 
a quasi-spherical, X-ray hot plasma. 

Another independent diagnostic of the accretion geom- 
etry is produced as the X-rays illuminate the disc. An iron 
fluorescence line and reflected continuum are formed wher- 
ever hard X-rays illuminate optically thick material (George 
& Fabian 1991; Matt, Perola & Piro 1991), but around a 
black hole these spectral features should be strongly smeared 
by the combination of Doppler effects from the high orbital 
velocities and strong gravity (Fabian et al., 1989). Modelling 
of the detailed shape of this relativistic smearing directly 
constrains the distance of the cool reflecting material from 
the black hole, while the amount of reflection and line con- 
strains the geometry. Results from this approach from low 
state spectra of both transient and persistent Galactic Black 
Hole Binaries are consistent with the geometry described 
above (Ueda, Ebisawa & Done 1994; Ebisawa et al. 1996; 
Zycki, Done & Smith 1997, hereafter Paper I; Zycki, Done 
& Smith 1998, hereafter ZDS98; Done & Zycki 1998). 

The inferred high state geometry is rather different. 
Both broad band spectra and detailed modelling of the re- 
processed component give a consistent picture in which the 
optically thick disc extends all the way down to the last 
stable orbit (Poutanen et al. 1997; Gierlihski et al. 1997b; 
ZDS98; Cui et al. 1997). The transition to the high state 
then seems to involves a decrease of r tr and consequently 
an increase of the disc thermal emission at the expense of 
the power law component. Further evolution to the very high 
state presumably requires a redistribution of the energy gen- 
eration so that the relative contribution of a coronal source 
increases. 

There is no current model able to compellingly explain 
all this evolution. Perhaps the best to date is that of Esin, 
McClintock & Narayan (1997), who go some way towards 
constructing the global dynamics as a function of mass ac- 
cretion rate. At low mass accretion rates there are least 
two thermally and viscously stable solutions of the accretion 
flow. One is the familiar optically thick, geometrically thin, 
SS disc, where the energy released is efficiently radiated via 



blackbody cooling. Another is a hot, optically thin, geomet- 
rically thick flow in which the radiative cooling is rather inef- 
ficient so that radial energy transport (advection) is impor- 
tant (advection dominated accretion flows, hereafter ADAF; 
see e.g. the review by Narayan 1997). For the hot solution, 
an increasing mass accretion rate, rh means that the flow 
density and, consequently, its radiative efficiency increases. 
Eventually, at some critical accretion rate rh = m cr it, the 
cooling becomes efficient enough to collapse the hot solu- 
tion down into an SS disc. Esin et al. (1997) calculate that 
this transition occurs for luminosities ~ 0.08Z/Edd, so they 
identify the low-high state change with the infalling mate- 
rial switching from an ADAF to SS disc. 

At least three of currently known transient black 
hole systems did not show this canonical behaviour. 
GS 2023+338 (V404 Cyg), GRO J0422+32 (Nova Persei 
1992) and GRS 1716-249 (Nova Oph 1993) remained domi- 
nated by the non-thermal component throughout the entire 
decline phase. For GRO J0422+32 and GRS 1716-29 this is 
plausibly because the system luminosities were never higher 
than ~ 0.1Z/Edd so that it never went above the low state 
(Nowak 1995; Esin et al. 1998, Revnivtsev et al. 1998), but 
this is clearly not the case for GS 2023+338, where the lumi- 
nosity reached almost the Eddington limit (Tanaka & Lewin 
1995). This system is especially interesting since it was used 
as an example for the ADAF-based models of SXT's in quies- 
cence (Narayan, Barret & McClintock 1997) , yet its outburst 
behaviour is apparently very different from that predicted 
by the extension of this model to the overall evolution of 
SXT (Esin et al. 1997). 

In this paper we present results of detailed spectral 
analysis of data obtained during the 1989 outburst of 
GS 2023+338 (Kitamoto et al. 1989), and use the relativis- 
tically smeared reprocessed component to constrain the ac- 
cretion geometry. The outburst was well covered by Ginga 
from its first detection by All Sky Monitor on May 22 un- 
til November 1989 (Tanaka & Lewin 1995) . In this paper we 
concentrate on the decline phase of the outburst after 1 June 
1989, leaving detailed discussion of the most dramatic "day 
in the life" of GS 2023+338 (30 May) to a forthcoming pa- 
per. For our analysis we have selected several data sets when 
the effect of photo-electric absorption (another unusual fea- 
ture of the source) was small and we aim at describing the 
evolution of the accretion flow geometry in the vicinity of 
the central black hole. We then compare this evolution to 
more typical cases of STX e.g. Nova Muscae 1991, and with 
the predictions of the Esin et al. (1997) model. 



2 DATA REDUCTION AND BACKGROUND 
SUBTRACTION 

The data were extracted from the original First Reduction 
Files in the usual manner using the Ginga reduction soft- 
ware at Leicester University. However, unlike most X-ray 
sources, the spectrum of GS 2023+338 contributes a signifi- 
cant number of counts above 24 keV. This count rate is used 
as a background monitor (Surplus above Upper Discrimina- 
tor, or SUD), but here it is contaminated by the source. The 
SUD rate minus the contribution from the source was esti- 
mated from off-source data taken before and after the source 
observation. We are able to do this because the satellite po- 
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Figure 1. Ginga light curve (1-20 keV) of GS 2023+338 since the first pointed observations on 23 May 1989 until November 1989. The 
count rate was corrected for background, deadtime and aspect. Dramatic flux and spectral variability was observed on 30 May. Flux 
variability occurred during each observation but the degree of spectral variability was changing with time. The upper envelope of the 
light curve (the 23. and 30. May points excepted) follows an approximately exponential decline with e-folding time ~ 40 days. 



sition in orbit, and thus the intrinsic background, depends 
on two parameters (excluding the phase of the 37-day orbital 
period of Ginga): the angle from ascending node of the or- 
bit and the longitude of the ascending node (K. Hayashida, 
private communication). We mapped the SUD rate from 
the uncontaminated observations against both orbital pa- 
rameters (dividing the data into 20 degree bins). For each 
data point from the observation of the source we calculated 
the orbital parameters and replaced the contaminated SUD 
value with that from the corresponding orbital parameters 
on the uncontaminated grid. We can then perform either 
a "local" or "universal" background subtraction as normal 
(Hayashida et al. 1989). We chose the universal background 
since there were no background observations suitable for the 
local method. 

To test the method we used the off-source (uncontam- 
inated) data using both the original and modified SUD val- 
ues. The r.m.s. deviations of the resulting spectra were less 
than 5-7 per cent over the 2-20 keV range. Secondly, we per- 
form a similar analysis on MCG- 2-58-22 (Nov 1989 obser- 
vation), comparing the "universal" background subtraction 
before and after inserting modified SUD values. The results 
from the real and modified SUD data were in agreement 
within the error bars, showing that the method is reliable. 



3 DATA SELECTION FOR SPECTRAL 
ANALYSIS 

The variability characteristics of these data were examined 
in detail by Oosterbroek et al. (1996, 1997). Here we show 
the light-curves (Fig. [lj) and colour-colour diagram (Fig. ^ 
for completeness. There is remarkable variability both in 
flux and in spectral shape (see also Tanaka & Lewin 1995). 
Oosterbroek et al. (1997) showed that the well defined track 
(marked Nn on Fig. ||) is consistent with being due to in- 
creasing photo-electric absorption. The left hand side of this 
track seems to end at a well defined point, plausibly because 
this represents the unabsorbed primary spectrum. We ex- 
tracted datasets from this point in order to try to constrain 
the underlying spectral shape since the beginning of June to 
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Figure 2. The colour— colour diagram of GS 2023+338 using the 
decline phase data i.e. after 1 June 1989. The well defined (al- 
most) horizontal track is due to photo-electric absorption with 
JVh increasing to the right. The turnoff point corresponds to 
JVjj ~ 4 X 10 23 cm -2 , when the spectrum above 5 keV begins to 
be affected. The circle marks a region where the June - July data 
for detailed spectral analysis were extracted from. The November 
data are those marked with filled triangles (see Table [lj) . 



the last Ginga observation in 1989 (November), as detailed 
in Table |l[ 



4 MODELS 

We use a variety of models for the data analysis. The soft 
component is described by a blackbody or a disc blackbody 
(i.e. multi-temperature composition of blackbody spectra; 
Frank, King & Raine 1992; Mitsuda et al. 1984) . For the hard 
component we use a power law, and its reflected continuum 
is modelled using the angle-dependent Green's functions of 
Magdziarz & Zdziarski (1995) as implemented in the pexriv 
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Table 1. Observations log for analysis of un-absorbed, decline phase spectra. 



label 


data set 


start date/time 


end date/time 


start date— 
30 May (d) 


exposure (s) 


1. 


June 3. 


89/154 3/06 16:36:31 


3/06 16:43:27 


4 


416 


2. 


June 4. 


89/155 4/06 12:11:59 


4/06 12:20:07 


5 


488 


3. 


June 10. 


89/161 10/06 01:58:00 


10/06 02:15:04 


11 


1024 


1. 


June 19. 


89/170 19/06 06:29:40 


19/06 07:46:40 


20 


1900 


5. 


June 20. 


89/171 20/06 14:05:00 


21/06 04:00:00 


21 


12180 


6. 


July 6. 


89/187 6/07 09:47:00 


6/07 10:10:28 


37 


1408 


7. 


July 20. 


89/201 20/07 01:00:00 


22/07 05:13:00 


51 


7650 


8. 


November 


89/305 1/11 00:20:00 


1/11 22:00:00 


155 


19064 



Table 2. Results of model fitting to the sequence of decline phase spectra. 



label 




fcT soft (keV) 


An (Rg) 


r 


N 




5 


R ln (R g ) 


X 2 /dof 


1. 
2. 
3. 
4. 
5. 
6. 
7. 
8. 


1.73 ±0.22 

1 54+ 026 
liO ^-0.20 
, 7 c+0.22 
L - ,o -0.24 

2 04+ ' 22 
Z ' U4l -0.26 

1 OO+0.30 
l - 2ii -0.34 

2 05+ ' 22 
Z ' UO -0.33 
1 1O+0-36 
J - x -0.25 
3.3 ±0.6 


n 99 +0.06 
U ' Z -0.09 
24+ 07 
n 97 +0.03 
u ' z '-0.04 

0.295 ±0.010 
37+ 12 

u -°' -0.04 

23+ 04 
u,zo -0.07 

43+ 22 
u, ^°-0.07 

30+ 04 


ll + °° 
7+82 
'-4 
o+5 
8 -3 

5 4+ 12 

1 6+ 10 
1 '°-1.3 

9.5+f 
6+ 07 

1+ 8 
'-—OA 


1.67 ±0.03 
1.70 ±0.03 
1.76 ±0.03 
1.72 ±0.03 
1.67 ±0.03 
1.64 ±0.03 
1.70 ±0.03 
1.93 ±0.05 


1 42+ 09 
^ -0.07 
i 91 +0.07 
1 - zl -0.06 
9 ot-+0.10 
z - od -0.13 

1.16 ±0.06 

0.67 ±0.04 
0.39 ±0.03 
(8.4±°;|) x io- 2 


0.70+0;- 
0.83«;» 

0.88+^ 
66+° 11 

u ' oo -0.07 

0.50 ± 0.08 

0.47 ±0.09 

0.37 ±0.07 

36+° 16 
U " 3D -0.10 


1+ 2 
1+ 2 
+o.i 

0+ 1 
1+ 8 

1 -0.9 

It? 
19 +30 
-8 

60 +<) " 
DU -45 


20+f 

30+fl 

17+1 

18+ 7 3 

36+?? 

sol- 

60+^ 

6001^65 


18.6/23 
21.6/23 
21.3/23 
23.1/23 
7.85/23 
22.8/23 
19.9/23 
26.3/22 



Nil - hydrogen column of uniform absorber in units of 10 cm 
D i n — inner radius of emitting disc as found from the amplitude of the soft thermal component, 
r — photon index of the primary power law. 
N - normalization of the primary power law (photons/cm 2 /s/keV) 
Q T - amplitude of the smeared, reprocessed component. 
£ - ionization parameter, Lyt/nr 2 . 
R ln — inner radius of the disc as found from the magnitude of the relativistic effects. 



model in XSPEC spectral analysis package (Arnaud 1996). 
This calculates the ionization state of the material by bal- 
ancing photo-electric ionization from the irradiating hard 
power law against radiative recombination, as in Done et 
al. (1992). The photo-electric absorption edge energies for 
iron were corrected from the rather approximate values given 
by Reilman & Manson (1979) to those of Kaastra & Mewe 
(1993). The recombination rate is a function of temperature, 
while the photo-ionization rate depends on the ionization 
parameter £ = Lx/nr 2 where Lx is the ionising luminos- 
ity, n is the density of the reflector and r is the distance of 
the reflector from the X-ray source. At a given temperature 
(which is chosen rather than calculated self-consistently in 
this code) then the ion populations are determined by £, 
with a rather weak dependence on the spectral shape (see 
Done et al. 1992). 

The iron Ka line corresponding to given ionization 
state, plasma temperature and spectral shape is computed 
self consistently in our code, rather than modelled as a sep- 
arate component. We use the Monte Carlo code of Zycki 
& Czerny (1994) to compute the line intensity and profile 
(which is broadened by Compton scattering in the reflector 
and a possible range of ionization stages of iron). The line 
is then added to the reflected continuum to give the total 
reprocessed spectrum. 

This total reprocessed spectrum can then be corrected 
for relativistic and kinematic effects which cause "smear- 
ing" of spectral features (e.g. Fabian et al. 1989). These ef- 



fects are calculated for a monochromatic line in the XSPEC 
diskline model so we use this, modified to include the light 
bending effect in Schwarzschild geometry, to calculate the 
spectral smearing, and convolve this with the total repro- 
cessed spectrum. The diskline model is parameterised by 
the inner and outer radii of the reflecting accretion disc, 
.Rin and R ou t, radial distribution of irradiation emissivity 
for which we assume Fi rr (r) oc and the inclination of the 
disc, i. In our basic models we assume ii out = 10 4 R g and 
(3 — —3 (see Appendix A) and fit i?j n and the amplitude 
of the reprocessed component expressed as the solid angle 
of the reprocessor as seen from the X-ray source normal- 
ized to 2-7T i.e. Q T = Q./2-k (so that Q r = 1 corresponds to 
an isotropic source located above a flat, infinite disc). This 
means that we do not assume any specific geometry which 
would result in a unique relation between Ri n , R ou t, P and 
f2 r . The models will be referred to as relrepr hereafter. 

We assume cosmic abundances of Morrison & McCam- 
mon (1983) with the possibility of a variable iron abundance 
in the reprocessor; we fix the inclination of the source at 
i = 56° and assume the mass of the black hole is 12 Mq 
(Shabaz et al. 1994). 

The spectral analysis was performed using XSPEC ver. 
10 (Arnaud 1996) into which all the non-standard models 
mentioned above were implemented as local models. 
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Figure 3. Unfolded spectra of the source (upper box in each panel) and x 2 residuals (lower boxes) for all data sets. Labels in upper right 
corners refer to given data set (Table g). Data points are plotted with (invisible) error bars. The solid histograms show the reprocessed 
component; the Fe Ka line and the Compton-reflected continuum are separated for plotting only. The dashed histograms are the primary 
power laws while the dotted histograms show the soft thermal component. Note that the flux scale in the last panel (8.) is shifted down 
by factor of 10. 



5 RESULTS 

The primary spectrum of the source is consistent with a sim- 
ple power law after 1. June (MJD-47000 = 677), so we can 
use this, together with the relrepr model and a soft compo- 
nent (disc blackbody) to describe the data. We assumed our 
basic reprocessing model, i.e. only one reprocessed compo- 
nent (cf. Section 3.1), with ionization and relativistic smear- 
ing allowed for. We fix R ou t at 10 4 R g and Fi rr (r) oc r~ 3 and 
we fit the amplitude fi r , the inner radius 7?i n and the ion- 
ization parameter f . Results of the spectral modelling are 
presented in Table and the spectra are plotted in Figure til 



The reprocessed component is significantly present in 
all data sets. It is very weakly ionized, and generally sig- 
nificantly smeared although not to the extent expected for 
a disc extending down to the last stable orbit at 6 7? g . We 
note that the two effects (ionization and relativistic smear- 
ing) can be separated given good quality data. The basic 
reason for this is that the ionization broadens the spectral 
features towards higher energy whilst the relativistic effects 
smear them around their rest-frame energy. Secondly, the 
strength of the Fe spectral features increases dramatically 
with ionization whilst relativistic effects only redistribute 
the photons. Figure H illustrates the above discussion show- 
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Figure 4. This plot demonstrates that the effects of relativis- 
tic smearing and broadening due to multiple ionization can 
be separated in our modelling. It shows contours of A% 2 = 
2.3, 4.6 and 9.2 as function of the inner disc radius R ln and the 
ionization parameter £(r = R; n ) (assuming £(r) oc r -4 ) for the 
June 19. data set. The best fit has low ionization (£ ~ 0.1) and 
R ln ~ 20, and there is no correlation between the two parameters. 



ing that there is no correlation between £ and R ln , even if a 
radial distribution of ionization is assumed. 

Figure |B^, shows how the amplitude of the reflected com- 
ponent decreases during the decline phase. Initially, the am- 
plitude of reflection is close to unity i.e. the cold, optically 
thick disc covers ~ 2ir solid angle as seen from the X-ray 
source, but this drops significantly as the source luminosity 
fades. The inner disc radius derived from the amount of rela- 
tivistic smearing is always incompatible with the last stable 
orbit in a Schwarzschild metric of i? ms = 6i? g (although 
with the caveat that i?i n is correlated with the form of irra- 
diation emissivity). The inner disc radius is consistent with 
being ant i- correlated with the amount of reflection, although 
the change in the amount of smearing is only marginally sig- 
nificant. 

These results agree with and extend those obtained in 
Paper I, where only a subset of present data was analyzed. 
They confirm that optically thick accretion material (the 
putative accretion disc) was present in the system through- 
out the decline, even though the source luminosity is clearly 
dominated by the power law component, rather than by di- 
rect thermal emission from the disc. 



6 COMPARISON WITH OTHER SXT 

The trend we see in the data of a gradual decrease of fi r 
possibly correlated with increasing Ri n is similar to that 
previously found in another SXT source, Nova Muscae 1991 
(ZDS98, Fig. 2). However, there is a marked difference in 
the spectral evolution of GS 2023+338 compared to Nova 
Muscae. The spectral index of the primary spectrum of 
GS 2023+338 stayed roughly constant during the decline 
phase, except for the last data set (November) obtained 5 
months after the outburst, when the power law was signifi- 



cantly softer than during the first two months. Conversely in 
Nova Muscae we found a monotonic decrease of F (harden- 
ing of the spectrum) with time, correlated with decreasing 
fir and increasing Ri n (see Figure [| ZDS98). A similar cor- 
relation between f2 r and V is also seen from GX339-4 (Ueda 
et al. 1994, Figure |) 

The results from Nova Muscae and GX339-4 can be un- 
derstood in terms of a phenomenological scenario where a 
retreating disc (increasing i?j n ) gives a smaller solid angle 
for reflection and a diminishing supply of soft photons for 
Comptonization, so leading to a harder X-ray continuum 
spectrum (Esin et al. 1997; ZDS98). The lack of such a cor- 
relation in GS 2023+338 is then a clear problem, implying 
that the observed change in solid angle for the reflecting ma- 
terial does not give rise to a change in the ratio of soft seed 
photons available for Compton scattering i.e. that the ob- 
served reflecting material is not the dominant source of soft 
seed photons for the Comptonization which gives rise to the 
hard power law. It is rather hard to imagine how this can be 
the case. If there is another soft seed photon contribution, 
perhaps from stronger cyclo/synchrotron emission, then the 
expected behaviour would be that the hard X-ray source 
would have a softer spectral index. GS 2023+338 instead 
shows a rather hard spectral index at all times. Perhaps a 
better explanation is that there is an additional contribution 
to the reflected component from more distant, cold material, 
which does not then contribute much to the soft seed pho- 
tons impinging on the X-ray source. We will now investigate 
such additional reprocessed components. 

6.1 Additional sites of reprocessing 

Additional material around the source is clearly present, as 
evidenced by the rapid spectral variability in the June 3. 
- June 19. data (see Figure []) which can be attributed to 
photoelectric absorption (Oosterbroek et al. 1997). The re- 
processed spectrum coming from this (presumably distant) 
material would not be relativistically smeared and its pres- 
ence in the spectrum could mean smaller amplitude of the 
smeared, disc reflected component. However, it is also prob- 
ably not optically thick to electron scattering. The spec- 
tra extracted from the extreme right (most absorbed) end 
of the 'Nh track' in Figure | give Nn « 4 x 10 23 cm -2 
(see also Oosterbroek et al. 1997). The Compton reflec- 
tion hump from such reprocessor is therefore suppressed 
but the iron Ka line is not, since the optical depth for 
photo-electric absorption at the Fe K-edge energy is ~ 1. 
We have used our Monte Carlo code to simulate reprocess- 
ing by such material including the iron line, assuming the 
reprocessor is not ionized. We then created a table model 
for XSPEC (parameterised by F) and repeated our spectral 
fitting of the smeared reprocessing including also this ad- 
ditional non-smeared component. The best fit contribution 
of this component is rather small, with fi r , n on-smcarcd Si 
0.01, but with 90 per cent upper limits (A^ 2 = 2.7) of 
fi r , non-smeared = 0.19, 0.22 and 0.19, for June 3., June 4. 
and June 10. data sets, respectively. However, this does not 
substantially change or broaden the allowed ranges of Q r of 
the disc-reflected, smeared component. Figure |7J shows the 
model fitted to the June 4. data set, with the contribution 
from the non-smeared reprocessing at its 90 per cent upper 
limit, and it shows clearly how the lack of a fully developed 
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Figure 5. Contours of A\ 2 = 4.61 (90 per cent confidence limits for two parameters) as function of the amplitude of the smeared 
reprocessed component, £7 r , and the inner disc radius, Ri n , derived from the relativistic smearing. Labels indicate different data sets 
(Table H), from the earliest (June 3.) to the last observation (November). Panel a shows the contours when only the disc-reflected, 
smeared reprocessed component is present, while panel b shows the same contours when an additional, non-smeared reprocessing is 
allowed for in the model. The amplitude of that non-smeared component is constrained to be < 0.2. Contours for only the first four 
data sets are plotted in b. Presence of the non-smeared component broadens the acceptable range of parameters of the smeared, disc 
reflection. 
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Figure 6. Contours of A% 2 = 4.61 (90 per cent confidence limits for two parameters) as function of the amplitude of reflection, £7 r , and 
the photon spectral index, T, for the decline phase spectra of GS 2023+338 (panel a; see Table ^) and (panel b) for the soft and hard 
states of Nova Muscae 1991 (ZDS98) and GX 339-4 (Ueda et al. 1994). The time arrow refers to Nova Muscae only. 



Compton hump in the optically thin reflection component 
cannot contribute much to the observed spectral harden- 
ing beyond 10 keV, so as much optically thick reflection is 
needed as before. 

Clearly, the observed (optically thin) absorbing mate- 
rial is insufficient to produce an excess reflected component. 
We cannot however exclude a priori a possibility of an addi- 
tional, optically thick reprocessed component being present 
in the data. Such component could come from e.g. an outer 



rim of the accretion disc, since the SS solution predicts that 
the disc flares in its outer regions. 

Solving the SS equations for vertically integrated, sta- 
tionary, non-irradiated disc we obtain O riIKm -smeared ~ 
H(R — i?out)/^?out ^ 0.03 at .Rout = Rtidai — 8 x 10° R g 
(for a = 0.01-0.5, rh = 0.01-0.05). Solutions with explicit 
vertical structure computed by Meyer & Meyer-Hofmeister 
(1982) give H/R(R out ) ~ 0.1. An irradiated disc would be 
expected to be more puffed up, but solutions of the vertical 
structure of irradiated discs depend on details of treatment 
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Figure 7. Model spectrum containing an additional, optically 
thin reprocessed component as fit to the June 4. data set. The 
amplitude of that reprocessed component is plotted at its 90 
per cent upper limit (the best fit value is factor of 10 lower). 
Such a component, due to reprocessing by cold plasma with 
Nn ~ 4xl0 23 cm -2 i.e. r es ~ 0.3, could be present in the data un- 
til June 19. The evidence for presence of such plasma comes from 
rapid spectral variability which can be modelled as photo-electric 
absorption (see Figure fel). Suppression of the Compton reflection 
bump due to t cs < 1 is clearly seen in the model spectrum. Be- 
cause of that, the required amplitude of the usual, optically thick 
reprocessing (possibly relativistically smeared) is still high. 



of the X-ray absorption. Meyer & Meyer-Hofmeister (1982) 
assumed that all incident X-rays are absorbed and obtained 
significant thickening of a disc, up to H/R(R out ) ~ 0.3. 
On the other hand, observational estimates of the optical 
(disc reprocessed) emission in Low Mass X-Ray Binaries 
also seem to imply H/R ~ 0.2 (de Jong, van Paradijs & 
Augusteijn 1996). 

We therefore added a cold (£ = 0), non-smeared re- 
processed component (continuum with iron Ka line) to our 
basic model and we repeated fits to first four data sets. In 
view of the above discussion we set an upper limit of 0.2 
to its amplitude. We found that a substantial contribution 
from such a component is indeed allowed by the data. Con- 
sequently, the contribution of the smeared component can 
be smaller. The results in the il r -Ri n plane for the relativis- 
tic reprocessing now allow for almost constant R[ n and Q T 
(Figure |b). 

Any reprocessed component coming from the distant 
outer rim of the disc should be constant or only slowly vari- 
able. Given that the reflected spectrum is hard, then this 
should reduce the amount of variability seen at high ener- 
gies, and around the narrow iron line. Oosterbroek et al. 
(1997) show that there is no obvious effect at the iron line 
energy for any of the decline spectra (although there is a 
marked constant iron line component during the outburst). 
However, this does not rule out the amount of unsmeared re- 
processing required here. Figure |^ shows the June 10th data 
set, together with the maximum contribution from the non- 
smeared reflection. Even though the line is intrinsically nar- 



' — _ so'urcr: 




Figure 8. Narrow spectral features in model spectra (e.g. an 
emission lines) are broadened by detector response, so their pres- 
ence does not manifests as a channel-to-channel change of e.g. rms 
variability in the observed spectrum. The upper panel shows the 
count spectrum (i.e. as registered by Ginga LAC) of the source on 
June 10 (data set 3.) showing the source, background and possible 
contribution from the model non-smeared reprocessed component 
(at it 90 per cent upper limit, cf. Figure ^). Lower panel shows 
the observed rms variability (computed from non-bgd-subtracted 
data) . 



row, it is broadened by the response of the detector. At the 
iron line energy the narrow component contributes only 3 
per cent of the counts, so its effect on the variability is small. 
Figure |^ also shows the observed 0.01-8 Hz rms variability 
spectrum of the June 10 data, showing that the variability 
pattern cannot rule out possible presence of an unsmeared, 
constant reflected component with fi rj , 10 n-smcaicti ~ 0.2. 

We are left then with the far from satisfactory situation 
where GS 2023+338 shows a very different behaviour to that 
of other SXT. Either it has a remarkable outer disc, which 
subtends a surprisingly large solid angle, or the source is 
somehow able to maintain a constant hard X-ray spectral 
shape despite there being a copious and changing amount 
of soft photons from the nearby accretion disc. 



7 DISCUSSION 

The results from the amplitude of reflection and amount of 
relativistic smearing imply that the disc does not extend 
down to the last stable orbit at 6 7? g . An independent esti- 
mate of the inner disc radius could in principle be obtained 
from the amplitude of the thermal emission, which is pro- 
portional to the emitting area. The amplitude is not always 
well constrained for our data because only the Wien cutoff is 
within the Ginga LAC band-pass for the observed tempera- 
tures. In those cases where the amplitude is well constrained 
the inner radius can be computed from 



A„ = D w 



K 

cos i 



Jcol 



[km], 
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where K is the normalization of the diskbb model in 
XSPEC, Dio is the distance to the source in units of 10 
kpc, i is the inclination angle and T co i/T e ff is the ratio of 
colour to effective temperatures, usually in the range of 1.5- 
1.9 (e.g. Shimura & Takahara 1995). 

The inferred inner disc radius {D\ n in Table ^ - N.B. 
values shown are computed without the colour temperature 
correction) is compatible with that inferred from relativistic 
smearing until June 19th. D lrl is smaller than i?j n afterwards, 
although a presence of a disk blackbody component of am- 
plitude corresponding to D ln — R[ n cannot be excluded, if 
its temperature was below ~ 0.18 keV, again due to limited 
band-pass of the LAC detector. 

The fact that D[ n and Ri n are in agreement is not easily 
compatible with strong constraints on ionization of the re- 
processor. Throughout the entire decline phase the reflect- 
ing material had to be very weakly ionized, with typical 
upper limit on the mean ionization stage of iron <FeV. 
Even assuming a strong radial distribution of ionization, 
£(r) = £(.Rin) (r/Rin)~ 4 , and allowing for the additional 
non-smeared and cold reprocessing (Section 6.1), we obtain 
the highest allowed Fe ionic stage of FeXVI. This corre- 
sponds to the reprocessor temperature < (2 — 4) x 10 K w 
20 — 40 eV, largely independently of whether iron ioniza- 
tion is in LTE conditions (e.g. Rybicki & Lightman 1979) 
or photoionization dominates (Figure 2 in Zycki & Czerny 
1994), and significantly smaller than the temperature of the 
observed soft emission. 

The observed (weakly significant) increase of £ in July 
and November data is not consistent with a stationary SS 
disc, since the latter solution predicts a weaker m depen- 
dence of the density n(rh) oc m 11//2 ° than Lx(rn) which 
presumably is oc m, hence £(rh) oc m 9//2 ° should decrease 
as m decreases. Similarly, radial dependence of the density, 
oc r- 15 / 8 is slower than that of the irradiating flux, oc r . 
This may mean that the disc is not in a steady state; in- 
stead the disc may have returned to cold, low viscosity state 
in its outer region where a mass build-up has begun, cut- 
ting off the supply of mass to the inner disc and causing 
departure from stationary solution for the inner disc. How- 
ever, solving the vertically averaged SS disc equations (for 
an un-irradiated disc) we find that Hydrogen recombines 
only beyond r ~ 10 5 R g when m ~ 10 -3 A#Edd, and so the 
stationary disc solution should be valid at the much closer 
radii where the reflection comes from. 

The power law spectral index in the last data set 
(November, 155 days after outburst) is significantly steeper 
than earlier (r ~ 1.9 vs. F ~ 1.7; Table 0). Since the source 
luminosity at that time was ~ 3 x 10~ I/Edd it may seem 
the source was close to quiescence and consequently its emis- 
sion was dominated by (Comptonized) synchrotron radia- 
tion, as in ADAF-based model of Narayan, Barret & Mc- 
Clintock (1997), and comparable to what was observed by 
ASCA. However, we do see the reprocessed component with 
Sl r ~ 0.3 so, consequently, the emission should be dominated 
by comptonization of the soft photons due to thermalized ir- 
radiation, if the reflection was indeed from the inner disc. Al- 
ternatively, again a significant contribution from outer disc 
is a possibility, with no feedback of soft photons to the X- 
ray source. However, the required magnitude of disc flaring 
is again rather larger than expected. We also note that if 
the outer disc contribution to reprocessing is dominant at 



the end of the decline phase, then the overall decrease of fi r 
has to be attributed to decrease of the smeared component, 
again rising the problem of the unchanging spectral shape. 



8 CONCLUSIONS 

Our results support growing evidence that the geometry of 
accretion flow evolves as the accretion rate changes during 
the decline of SXT. The evolution involves a retreat of the 
inner, optically thick disc as indicated by both decreasing 
amplitude of the reprocessed component and the relativistic 
smearing becoming less or insignificant (Paper I; ZDS98). 
These variations of geometry, observed also in persistent 
source Cyg X-l (Done & Zycki 1998) are usually correlated 
with spectral states: small inner disc radius corresponds to 
soft spectra whilst large radius corresponds to hard spec- 
tra (Nova Muscae 1991, ZDS98). This correlation is easily 
understood in terms of a scenario where the availability of 
soft seed photons for Comptonization, determined by the 
inner disc radius, regulates the Compton cooling and, con- 
sequently, the spectral shape (Esin et al. 1997). 

However, while GS 2023+338 showed the same reced- 
ing inner disc, its hard X-ray continuum spectral shape re- 
mained largely constant. This poses serious problems to any 
model involving Comptonization as a source of hard X-rays, 
and in which changes in the spectral properties are linked 
to the geometry of accretion. 
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APPENDIX A: RADIAL DEPENDENCE OF 
IRRADIATION EMISSIVITY 

The irradiation emissivity follows closely the r -3 depen- 
dence if the reflecting disc is outside a central source of 
radiation, except very close to the source. Consider a spher- 
ical source, centered on the black hole, of maximum radius 
Ro, with the disc extending outside it, i.e. from r = Ro to 
Rout 2> -Ro- We allow for a flaring disc with height described 
by an arbitrary function h(r). Introducing spherical coordi- 
nate system in the usual manner (see Figure Al) we have: 



S = [p sin 6 cos <f>, p sin 6 sin <j>, p cos 6 
D = [0,r,A(r)]. 

Vector normal to the disc surface is 

-± 1 

dr 

Distance between S and D is the length of the vector from 
S to D, I = \x\, 

I 2 = p 2 + r 2 - 2pr sin 6 sin cf> + h(r) [h(r) - 2p cos 6], ( Al) 
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Figure Al. Geometry of disc irradiation. The X-ray source is assumed spherical, with outer radius Ro- The disc extends from Ro to 
Rout 2> Ro- S is an arbitrary point within the source while D is a point on the disc surface. Note that the X— ray source is plotted in 
perspective but a cut in y-z plane for the disc is used. 



while the cosine of the illumination angle is 
COS7 = n ■ a;/(|n||a;|) = 



4^ (r — p sin 9 sin < 



+ p cos 9 — h(r) 



(A2) 



I yjl + (dh/dr) 2 

Contribution to the irradiation emissivity is thus 

6F ilr = dV cos j /I 2 , (A3) 

where dV = p 2 sin 9 dp d9 dcj> is the volume element of the 
source. Neglecting radiative transfer and relativistic effects 
within the source, we obtain the emissivity by integrating 
SFirr over the source, 

Firr(r) = J f(r)SF irl , 

where f(r) describes the distribution of volume emissivity 
within the source. 



Examples of Fi rr (r) are plotted in Figure A2. As can be 



seen close to the source the emissivity can deviate from the 
r~ 3 dependence. It is generally steeper (r - ^ with /3 > 3) due 
to the finite size of the source. The magnitude of deviation 
depends on luminosity distribution within the source. The 
stronger the central concentration of luminosity, the weaker 
the deviations. 

A special case of a point source located on the symmetry 
axis, above a disc is recovered from (Al) and (A2) assuming 
polar angle 9 — 0. Assuming further a flat disc (h(r) = 0) 
we obtain well known formula, 



Fi rr (r) (X 



1 



(r 2 + h 2 f /2 



for r > h s 



where h B is the height of the source above the disc plane. 

Another frequently quoted formula for the irradiation 
temperature in a flaring disc (SS, van Paradijs 1996, King 
et al. 1997), results from our formulae (Al) and (A2) in the 
limiting case of p = 0, putting I — r (i.e. it is equivalent to 
assuming r ^ Ro) and neglecting dh/dr in denominator of 
(A2). 
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Figure A2. Irradiation emissivity, r s Fi IT (r), for a central, spher- 
ical source illuminating an external, flat disc. Inner source ra- 
dius is assumed 6R S , outer one 20i? g , from where the disc ex- 
tends outwards. The emissivity is somewhat steeper than r~ 3 for 
r < 100-Rg and its shape depends on how strongly the luminos- 
ity inside the source is concentrated towards the centre. The solid 
curve shows the emissivity for a uniform source, whilst the dashed 
curve is for the luminosity per unit volume r~ 4 (1 — 6/r). In 
the 20 — 40-Rg range the approximate power law fit to Fi tt (r) has 
slope equal fi = 3.5 in the former case, while = 3.3 in the latter. 



Our formulae (A1)-(A3) could also in principle be ap- 
plied to a disc extending inside the source. A simpler way 
to find the irradiation emissivity in this case is to notice 
that Fi TI (r) follows the local electron emission, vertically 
integrated over the hot source, q~ . This generally follows 
the viscous energy generation in the hot source which how- 
ever is not necessarily described by the usual dependence 
q + oc r~ 3 . First, the ratio of q^ ot /qt t can be a function of 
radius (Witt, Czerny & Zycki 1997). Second, a fraction / of 
9hot can ^ e advected, and / itself can be a function of radius. 
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Third, Compton cooling by non-local soft photons can have 
a strong effect (Esin f997), especially in the transition re- 
gion. Since the non-local cooling enhances electron emission 
(Esin f997), we can expect Fi Tr (r) to be steeper than r~ 3 
i.e. p > 3 close to the inner disc radius. All the above effects 
are model-dependent and need to be properly modelled. 
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